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^L(. Abstract 

^ ■ 
^Q ' The motion of a classical spinning test particle in the field of a weak plane gravitational wave 

is studied. It is found that the characteristic dimensions of the particle's orbit is sensitive to the 

p^ ' ratio of the spin to the mass of the particle. The results are compared with the corresponding 

motion of a particle without spin. 
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O ■ 1 Introduction 

O 

^^ I One of the problems of great importance in contemporary physics is the detection of gravitational 

pis waves which would profoundly affect both physics and observational astronomy [Q]. Like other kind 

I ■ of waves, one may think of gravitational waves as wiggles propagating in some background. In fact 

{^X)! these waves are wiggles of curvature propagating in a fixed background with the speed of light. The 

corresponding space-time metric may be written as g^u = g^ + h^y in which g^ stands for the 



background metric and h^^, represents the gravitational wave field. In practical situations where the 
detection of these waves is our prime concern, one needs to consider only weak, plane, linear, waves 
C^ I in an essentially flat background because the amplitude of the waves crossing the solar system are 

very small even though the energy flux may be enormous [y]. We are thus led to the decomposition 
of the space-time metric as g^^ = rj^y + /i^,y where rj^y = (—,+,+,+) is the usual Minkowski metric 
and \h^y\ ^ 1. It then follows from general relativity that h^y satisfies a generalized wave equation 
in three dimensions and there are only two components for h^y corresponding to two independent 
polarizations of gravitational waves usually denoted by /i+ and h^ [§1- Polarization of these waves 
contains valuable information about the sources generating them [^. 

Gravitational waves alter the curvature of the region of the space-time they cross and hence affect 
the motion of particles. This problem has been studied extensively in the past e.g. see [^-|§. The 
basic idea is that free particles move on the geodesies of the basic field g^y. The physical effect of a 
passing gravitational wave is that particles acquire a relative acceleration during the passage of the 
wave and hence gain a relative velocity after the wave has passed. The measurement of this relative 
velocity would hopefully provide a way to detect gravitational waves. 

Grishchuck Q studied the problem of the motion of a particle in the field of weak plane grav- 
itational waves with different polarizations and calculated the trajectory of the particle as viewed 
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from a locally inertial frame, that is, from the reference frame attached to a nearby particle. These 
trajectories, when projected onto the plane perpendicular to the direction of the wave propagation, 
become ellipsis or segments of straight lines whose dimension and orientation depend on the initial 
position of the particle and the polarization of the wave. By a suitable combination of the particle's 
initial position and wavelength of the wave, it is possible to find invariant dimensionless parameters 
which characterize the dimension of the trajectories. 

In this paper we consider the problem of the determination of the motion of a classical free 
spinning particle in the field of a weak, plane gravitational wave with arbitrary polarizations. Here, 
several interesting phenomena will appear. The particle's spin couples to gravity and this makes the 
spin dynamical. Also, the trajectory of the particle deviates from geodesies of the field 5^,^ as a result 
of a Lorentz-like force which originates from the spin-gravity coupling. This problem was previously 
considered from a different point of view in 0] where a supplementary equation was proposed causing 
the spinning particle to move along a geodesic. 

In what follows we first review the equations of motion of a spinning particle in a general space- 
time, sometimes known as Dixon's equations [^], and derive their weak field approximation. We then 
apply them to the case of a spinning test particle in the field of a gravitational wave and compare 
the results with those of a test particle without spin. 

2 Dixon's equations 

The motion of an extended body in a given background may be studied by treating the body via 
a multipole expansion. This expansion is based on the following procedure. The body sweeps out 
a narrow tube in space-time as it moves. Let L be a line inside this tube representing the motion 
of the body. Denote the coordinates of the points on this line by X" and define 5x" = X" — x", 
with x" being the coordinates of the points of the body. Let T^^ be the energy-momentum tensor 
describing the body and consider the quantities / T^^dV , J Sx^T'^'^dV, J dx^Sx^'T^^'^dV, ■ ■ ■ where the 
integrations are carried out on the three-dimensional hypersurfaces of fixed time X^ = const.. These 
are in fact the successive terms in the above mentioned multipole expansion. If we only keep the first 
term, that is, if we assume that all the integrals containing dx^ vanish, this will describe a single-pole 
particle. Similarly, we may consider a pole-dipole particle for which all the integrals with more than 
one factor of 5x^ vanish. Higher order approximations may be defined in a similar way. In the limit 
of ||(5x"||^0 with ||.{| denoting the norm with respect to the induced metric on the three-dimensional 
hypersurface of constant time, the approximation becomes exact. Thus a single-pole particle is in fact 
a test particle without any internal structure. A pole-dipole particle is a test particle whose internal 
structure is expressed by its spin, an antisymmetric second rank tensor defined by 

dx^'T^" - (5x^T°'^') dV. 



The equations of motion are then obtained by applying Einstein's field equations together with the 
conservation of the energy-momentum tensor, V pT^^ = 0, describing the body. For a single-pole 
this leads to a free particle moving along the geodesies associated with g^y. This is the same as that 
obtained by Einstein and co-workers, who only made use of the field equations without resorting to 
the dynamical equation. Using the above procedure, Papapetrou [^ derived two sets of equations for 
the motion of a pole-dipole particle and its spin 
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'' = -W.Xo^''^^' (1) 



9 ii\p 



s'^^ = pt^v" - fv^ (2) 

in which v^ = — ^^^ is the four- velocity of the particle along L parameterized by r with v ■ v = —1, 
R^^Xo ^^ ^^^ curvature tensor, p^ = J T^'^dV represents the particle's four-momentum along L and 



.,. ^ ^ + r",„„V 



ds^"'iT) 



+ T^,v^sP'' + T\ys^'P 



are covariant derivatives. Equation (||) expresses the coupling of the spin with curvature, and for 
either a flat space-time or a spinless particle reduces to the geodesic motion p^ = 0. Also, equation 
(|2[) is in fact an equation showing the balance of the particle and the field angular momenta. Unlike 
special relativity, p'^ and v'^ are not generally proportional to each other. Equations (|l]) and (|2|) 
alone are not adequate to give a solution and three of the components of s^'^ remain undetermined 
by these equations. This may be regarded as a result of L not being yet fixed. To remedy this, 
several supplementary conditions have been proposed. A widely used condition proposed by Pirani 
1^] is v^s^'^ = 0. However, it was shown in the second reference in [|^ that this condition leads to 
unphysical solutions where the particle can admit helical paths in flat space-times. Also, it does not 
uniquely specify a world line through the body [|14| ]. Here, we adopt the following condition |jl^ 

PmS^' = 0. (3) 

This choice does not suffer from the appearance of unphysical solutions mentioned above and has 
other advantages discussed in |0, [T2| and p^ . Equations ([l|) and (^) were derived by Dixon |11, 12] 



in a covariant manner and equation (|3D was adopted as the supplementary condition. These equations 
have also been studied from different angles in |]l^ . The post-Newtonian approximation of equations 
(|l]) and (0) together with v^s^'^ = as the supplementary equation was studied in |15]. From Dixon's 
equations we deduce the following relations 

p^p^ = const. = —m , TiS^" s^u = const. = s . 

Alternatively, equation (P) enables us to express the particle's spin by the four-vector 

2m 
where ■q^'^ = ^^g/'^^^K ^g ^j^^ alternating tensor with gO^^S _ ^ ^^^ g ^g ^j^g determinant of the 



metric. It then follows that Pfj,s^ = and s^s^ = s^. Under certain conditions [16|, Dixon's equations 
may be thought of as the classical limit of the Dirac equation in curved space-time. 

To find the particle's trajectory it is necessary to first find its four-velocity. Since there is no 
equations of motion for this purpose, one has to find the velocity from Dixon's equations in an 
indirect way. One such relation is |17| 



{p ■ p)v>' = {p-v) (p^^ - _3^E]^^pl£!Zl!l\ . (4) 

Unlike for p^^ there is no guarantee from the equations of motion that v^ remains time-like. In fact 



Tod et.al. |17] have solved Dixon's equations for the case of the motion of a spinning test particle 



in the equatorial plane of a Kerr black-hole with the particle's spin perpendicular to this plane and 
found that there are orbits for which v^ becomes null or space-like. 

As Dixon's equations are reparametrization invariant, we have to fix a gauge before embarking 
on finding solutions for them. It is convenient to choose 

1 

— p ■ V = —1. 
m 

This is the gauge in which the instantaneous zero-momentum and zero-velocity frames are simulta- 



neous 



3 Spinning particles in a weak field 

It is generally difficult to solve equations ([l|)-(0) analytically and to our knowledge, the only solution 
that has been found up to now is the previously mentioned solution by Tod et.al. However, for the 
case of weak fields in a fiat background where the metric is of the form g^^i, = r/^,y -|- h^i,, it is natural to 
assume that the change in the particle's velocity and spin are of the same order as hfj,^,. Therefore, we 
may decompose the particle's four-momentum, velocity and spin tensor as p^ = Pq+tt^, v^ = Vq+w^^ 
and s^"^ = Sq'^ + a^^ with vr'^, w^ and a^ being quantities of the same order as h^y and pg, Vq and 
Sq^ are the initial (before the passage of the wave) four-momentum, four-velocity and spin tensor of 
the particle respectively. The particle's trajectory can then be found from — -^^^ = v^ subject to 
suitable initial conditions. Similarly we may express the particle's spin vector as s^ = Sq + a^ . In 
terms of these decompositions, the zeroth order equations become 

Equation (|5|) is trivial: in the absence of any field the particle's four momentum and spin tensor 
are conserved. Equation (^) eliminates three of the components of s^'^ and equation (^) may be 
used to express Sq^ in terms of Sq. For example, for p^ = (m, 0,0,0) these lead to Sq^ = and 

For the metric 5^,^ = g^ + h^y , the curvature tensor and the connection may be written as 

n _ 7->{0) I 1 Y^i _ -p(0)At I M 

^^Uk\ — ^f_iUKX ' ' f^UKX ailU i ,yf^ — i^ UK ' luK- 

However, the background quantities R\1^x and F l^ vanish for the flat background and the flrst order 
form of equations (|I|) , (Q) and (^) reads 

da^'^ir) _ 1 ^ „u ry ^>i.p ^.k\ ^,Ai ,,\,.pv 

JqV 

1 .., . 1 



dr 2m 



rp.nxv',v^,s^,'sl^ - i^pV^s^S" - /^--^ (9) 



^'{r) = -vr^(r) + ^r.Ap.^^o'«^«r (10) 



with the supplementary equation written as 

r/,,(.^o''^vr^(r) + p^a^'' {r)) + h^.p^,^ = 0. (11) 

The corresponding expression for a^ is 

<^'{r) = ±e^^,^ipPa''' + t^^s^,' + 2p^,So\h\ + p^.s'^'h'' ,). (12) 

Therefore, given an initial data set Vq = {vq,Vq,Vq,Vq) and Sg = (sj], sj, Sq, Sq), the orbit of the 
particle is given by 

X^{T)=X^{r)+x^{r) (13) 

and its spin components are determined by solving equations (|8|)-(p|) subject to suitable initial 
conditions. Here, — ^^^ = Vq and '^J^ = w^^. 

4 Gravitational waves 

In this section we apply the procedure described in the previous section to a spinning particle in the 
field of a weak plane gravitational wave propagating in the z-direction in some coordinate system. 
For such a field, the metric may be written as |p 

ds'^ = -df + (1 + a{t - z))dx^ + (1 - a{t - z))dy'^ + 2b{t - z)dxdy + dz^ (14) 

with a{t — z) = h^ sin(u;(t — z)) and b{t — z) = hx cos{uj{t — z)). A wave with linear polarization is 
described by /i+ = or /ix =0 and a circularly polarized wave is characterized by /i+ = ±/ix . 

For this metric, it is straightforward but lengthy to solve the relevant equations. The resulting 
expressions for the trajectory of the particle and the evolution of its spin are too complicated to be 
included here, but are given for x^i"^) iii the Appendix. However, they become manageable for the 
physically interesting case of the motion of a particle initially at rest with respect to the coordinate 
system in which the metric (14) is expressed. In this frame, Vq = 1,Vq = Vq = Vq = and an initial 



spin orientation may be taken as (0,So,Sg,So) . For such a particle the trajectory (with respect to 
the coordinate system {t,x,y,z) and the initial conditions x^i^) = w^{0) = cj^(O) = 0) is given by 

t{T)=T (15) 

x{t) = -- — ^sl{slujhx + mh+){uJT - sin(a;r)) + - — i^s\{s1ujh+ + m/ix)(l - cos(L<jr)) (16) 

1 1 

y(r) = -^^-^sl{sliohx + mh+){LOT - sin((jr)) - ■^^-^sl{slu}h+ + m/ix)(l - cos{u)t)) (17) 



1 1 

-^sl{slujhx + mh+){ujT - sin(a;r)) - — - 



z(t) = — 2SoSow/ix(u;t -sin(wT)) + -— 2u;/i+((so) - (sg) )(1 - cos(a;r)) (18) 

and the equations for spin evolution are 

s\r) = ^K {{slf - {slf) sin(^T) - -slslK{l - cos{u;t)) (19) 

s^ ('^) = •So - 2^lh+ sin(a;r) + -Sq/ix (1 - cos{u:t)) (20) 

s^{t) = sl + -slh+ sin(a;r) + -sj^x (1 - cos(a;r)) (21) 

s\t)=sI (22) 



We can use equations (2C)-(|2^ to study the evolution of the particle's spin. From these equations 
it is clear that the spin evolution does not depend on Sq, nor is this component affected by the wave. 
Therefore, the initial component of the spin along the direction of the wave propagation is unaffected 
by the wave. The components s^(t) and s'^{t) depend on both Sq and Sq giving rise to the spin 
vector precessing about the z— axis. The rotation of the spin vector is very small and as equations 
(pOJ-^) show, they are of the order of h^s or hxS, where s is the relevant component of the spin. 
The appearance of a nonzero Sq component is related to the fact that the particle is not at rest with 
respect to this coordinate system. 

The particle's trajectory with respect to the frame (t, x, y, z) can be determined from equations 
([l5|)-(|T8|). Putting Sq = Sg = and Sg = s in these equations result in 

(t(r),x(r),y(r),z(T))=(r, 0,0,0) (23) 

which means that the particle remains at rest at its initial position. This behaviour is the same as 
that of a particle without spin. This has no invariant geometrical meaning and may be interpreted 
as follows: the coordinate system {t,x,y,z) is the one always remaining attached to the particle. 
The physical effect of the wave is on the relative motion of the particles. To study this we switch to 
(i,x,y,z), a locally inertial frame associated with the worldline x(r) = y(r) = z{t) = given by [g] 

t = t--^{y'-x')--xy (24) 

(25) 

(26) 

(27) 

with a, b, da/dt, and db/dt being evaluated at x = y = z = 0. In this way, we are able to study the 
relative motion of two nearby particles, one with spin and the other without spin. 

The motion of a nearby particle without spin initially at the position {li,l2, 0) with respect to this 





Ida 2 2\ 9b 


t 


-Adt^y ") dt'^y 




1 1, Ida Idb 


X 


-\ — ax by t:-xz H r^yz 




2 2 ^ 2dt 2dr 




1 1, Ida Idb 


y 


ay bx T^yz -\ r^xz 

2 ^ 2 2dr 2dt 




Ida . 2\ f^^ 


z 


-Adt^y ") dt'^y 



frame was studied in [19|. The motion of a spinning particle with its spin parallel to the direction 
of the wave propagation with respect to this frame is the same as a particle without spin. Let us 
now consider the motion of a particle with its spin perpendicular to the z— axis. We thus set, say, 
Sq = s and Sg = Sg = with the initial position x = li,y = I2, z = 0. We then obtain from equations 
(|lD-(|D and (ID-dH) 

t(r) = T-'^{ll-ll)cos{uJT) + '^hl2sm{ujT) (28) 

x(t) = /i(l + -— sin(ujT)) ^l2 cos(ujt) + --^(1 - cos(a;r)) (29) 

2 2 2m 

y{T) = hil - — sin(u;T)) - -^h cos{ujt) - —^{ujt - sm{ujT)) (30) 

z{t) = — (^2 -ii)cos(cjr) -— /i/2sm(a;r) -— ( — ) (1 - cos(a;T)). (31) 

4 2 2 m 

This set of equations determines the trajectory of a spinning particle initially at rest at the point 
(^1, ^2; 0) with respect to a locally inertial frame associated with the origin. In general, this trajectory 
is in the form of a helix whose size is sensitive to the ratio of the particle's spin to its mass. In figure 
1 the trajectories of a spinning particle and that of a particle without spin are shown for ^ ~ 10"^ 



m. We note that in a system of units in which c = 1, this ratio for an electron is of the order of 
-|^Q-i3 j^ g^j^j £qj. g^j-^ object hke the Earth is 10^^ m. For — ~a;~^~A the spin contribution to the 
displacement in the z— direction is of the same order as those for the x and y— directions. Also, it 
is worth mentioning that for — ~/i~/2i the contribution of the terms involving the spin in equations 
(|29|)-(|3l|) are as important as the terms representing the motion of a spinless particle. Therefore, the 
ratio ^ seems to play a remarkable role in this solution. 

One may estimate the displacement of a test particle using equations (p8|)-(|3l|). During one period 
of the wave, i.e. ujt = 2tt , a particle initially at (li, h^O) would reach the space-time point 

i = T-\^h+{ll-ll) 

X = h--hxl2 

y = h - -f^xh - T^h+ — 
2 m 

z = --h+uj{ll-ll). 

corresponding to a displacement of the order of hi with / = Jl\ + Z^- Choosing the set of parameters 
u) ~ 10^ Hz, /i ~ ^2 ~ :^ ~ 10~^ m and hj^ ^ h^^ ^ 10~^^, the displacement would be of the order 
of 10~^^ m. The above values for oj and /i+, h^ stem from the expectation that for the gravitational 
waves accessible to earth based detectors a; ~ 1 — 10^ Hz with a typical amplitude of /i ~ 10~^^ [^. 
This value for uj together with ujt = 27r corresponds to r ~ 10~'^ sec which is also a typical duration 
for short bursts of gravitational waves Q . 

One would hope to be able to detect such gravitational waves by the future sensitive detectors 
such as LIGO [l| currently under construction. In general, these displacements are small but could 
be cumulative and hence likely to be observed. 

5 Conclusions 

We have studied the motion of spinning particles in the field of gravitational waves. The solutions 
we have found are free of the problems mentioned before, namely, the violation of the normalization 
condition v ■ v = —1. For a spinless particle, the velocity and hence its orbit depends on its initial 
position. For spinning particles, there is a contribution to the change in the trajectory from the 
coupling of the spin to the curvature. This is independent of the initial position and rather depends 
on the initial orientation of the spin. The amount of deviation from the geodesic path depends on 
the ratio ^. For a certain value of this ratio, the displacement caused by the spin-gravity coupling 
could be as important as that caused by the usual tidal force of gravity waves. Although the particle's 
displacement is very small, one may hope to have it amplified in the course of the wave passage, e.g., 
by a suitable arrangement of reflectors. 

In this paper we have only considered monochromatic waves i.e. those with a single oo. In 
principle it is possible to generalize this procedure to the case of gravitational waves which are not 
monochromatic [||]. Our approximate solutions may also provide some hints for finding exact solutions 
to Dixon's equations for a general plane wave. For example, we can investigate whether there are 
exact solutions to these equations corresponding to the motion in certain planes with specific spin 
orientation. Indeed, equation ([2^ ) is such an exact solution. This has subsequently been verified in 



[18 1 where a class of exact solutions along whose path the initial orientation of the spin of the particle 
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Figure 1: The path of a spinless particle (sohd hne) and a spinning particle (dashed line) in the 
field of a gravitational wave. They are originally located at the same position (-1,-1). The values 
of the parameters used here are lj~1 kHz, /i4_~/ix~10~^^, /i ~ ^2 ~ 10~^ m and — ~10^"' m. 
For convenience, we have adopted a system of coordinates given by X = 10^^ (x — 10~^) and Y = 
10^^(y — 10~^). The displacements in the Z direction are very small and thus are not shown in the 
figure. 



remains unchanged, is presented. It should also be interesting to study the cumulative effects of the 
waves on the motion of particles |19|. 

It is worth mentioning that within the context of the perturbative framework presented here, we 
have also used Tulczyjew's description (a special case of Papapetrou equations together with Pirani 
supplementary condition) |£^]. We have found that in some special circumstances, e.g. when the test 
particle is initially at rest, the spin evolution is exactly the same for various supplementary conditions 
mentioned above. However, in general the trajectories turn out to be only slightly different. The 
general case of Papapetrou-Pirani will be reported elsewhere |21|. 
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6 Appendix 



In this appendix, we give the components of the trajectory of a test particle [22| in the field of a 
gravitational wave specified by the metric (0). To make the expressions manageable, the following 
substitutions have been made: si, S2) ■S3,p, k stand for s\^s\.,s\.,{a — 6)iOT — sin((a — 5)ujt) and 
1 — cos((a — 6)ujt) respectively with {vq,Vq,Vq,Vq) = (a,/?, 7,(5). These components together with 



those of the spin (not presented here) satisfy equation (|5). Also, using equation (|1 
to equations ([l5|)-([l8|) for the set of initial parameters discussed in section 4. 



they reduce 



/(r) = {p [u;2/i^ (s?(a^/37 - 3a'^(3-fS - afi^-^ + a(3-/^ + 3ap-f6'^ + (3^-fS - l3-f^6 - (3-/6''^) 



,3^2 



+ siS2{2a^S-a-'P 
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+ 3a^p^-f - a^-i^ - 6aS(52 - 60/3^5 + 2a-i^5 + Aa^iS^ + 3/3^7(^2 _ ^3^2 _ ^^4^ 

+ s^(q3/37 - 3a2/?75 + a/3^7 - a/37^ + 3q/37(5^ - /J^*^ + h^^ - l^lS^) + S2S^{-a'^(5 

+ 4a3/3(5 - a2/?3 + 2,0^13-^"^ - Qa^ps^ + 2ap^S - Qa/3-f^6 + Aaf55^ - p^6^ + 3p-f^6^ 

- P6*) + sl{-2a^p-f + 6a^/37(5 - 6ap-fS^ + 2/37(5^)) + mh+uj{si{-a^-f + q7(5 + 2/3^7) 



+ S2{-a^P + a/3(^ + 2p-f^) + S3(-2a/37 + 2/37(5)) + 2m'^hy, a/37 
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x\t) = {p[uj^h,,{sl{a^P^-f-2a(3^-f5-(3S + Ph^ + PhS^) + siS2{-a*P + 3a^(36 
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a2/?3 - a2/?72 - 3a2/352 - a/3^5 - 0/3725 + a/35^ - 2/3^2 ^ 2/37^ + 2p-f^5^) 



2a/37^ - 4^/3752 - 2/3^75 + 2/37^5 + 2/3753 ) + s^(q^7 



+ siS3(2a2/375 + 2Q/3S 

- 0^75 + Q2/327 - 2a2-y3 _ Q,2^^2 _ Q,^2^j _^ ^^3j _^ ^^j3 _ ^2^3 _^ ^5 _^ ^3^2^ 

+ S2S3{-a^ + 2a*5 - a^p^ + 30^72 + 2a^p^6 - 4a^-f^6 - 2a^d^ + 2aP^-f^ - ap'^5^ 



,2„3 



2a7* - ^7252 + ad"" - 2p'Y6 + 27^5 + 2725^) + Sgl-a^ + 2a'^75 - a^P^j + a'j' 
2aP^-f6 - 207^5 - 2075^ - p^-f6^ + 7^52 + 75^)) + mh+uj{si{ap-f) + S2{-a^ + a^S 



"7^) + '53(-a^7 + a7'^)) + m.^^xa^7 +n h+u?{sl{-a^ P5 + a^P^i"^ + 2a^P5'^ 

ap^5 - aPj^6 - ap6^ - 2P^-/^ - p^d^) + siS2(-aS + "^5 + 3a2^2^ ^ ^2^3 ^ ^2^^2 

- a/3275 - 07^5 - 075^ - 4/^273 _ 2p^-f6'^) + sissia^ - 2a^b - o? p^ - o?-i^ 

+ 202725 + 2q253 + 40/3272 + 30/3252 - 07252 - o5^ - 4/52^2^ _ 2/j2^3) ^ sl{-a^P 

+ 2o?pb + 3o2/?72 - ^2/352 - 20/3725 - 2/37^ - /?7252) + s2S3(-4o3/37 + 6o2/?75 + 40/37^ 

- 4/37^5 - 2/3753) + s^(o^/3 - 2o?pb - 2o?p-i^ + Aap-i^b + 20/35^ - 2/37252 - /?5^)) 

+ ?Ti,/ixi^(si(o — o 5 — 0/3 ) + S2(— 0/37) + 53(0 P — aP6)) — m h^ a P > /{2m uoa ) 



X\r) 



-xHr 



-T,/3^7, si ^ S2) 



;^3(^) = |p [/ixL<j2(sf (-0^/37 + 3o3/375 + 02/3^7 + 02/373 - 302/3752 - 30/3375 - 0/3735 + 0/3753 



+ 2/33752) + siS2(2o6 - 60^5 - 3oV^ - 3oS^ + 60^52 + 7a^P^5 + 703725 - 2o353 
+ o2/34 + 2a^pW - ^c^pH^ + 0274 - 5027252 - 0/345 - 60/32725 + 0/3253 - 07^5 
+ 07253 + 4/327252) + siS3(oS - 4oS<5 - "3/327 - 0373 + 603752 + 602/3275 



.2„3 



+ 20^7^5 - 40^75^^ 



90/32752 - 07352 + 0754 + 4/52^53) + s^(-oV7 + 3o3/375 



- a^l3^ - a^p-f^ + 6a^f36^ + 2a^f3^6 + 602^72,5 - 40^/3(53 - a^'-^d^ - dap-f^d"^ + a/35^ 

+ 4/372(53) + s^(-2a3/375 + 6a'^p-f6^ - 6ap-/6^ + 2/37^^)) + m/i+L<j(si(-aS + a%5 + 2a/3S) 

+ S2(-a^/3 + 0^/3(5 + 20/37^) + S3(-2a2/J^ + 2ap-f6)) + 2m^h,^Q^(3\ + k \h+u:'^{s\{-a^ 

+ 2,a^5 + 20^/32 + aS^ - Sa^-J^ - ba^fi'^6 - 2a^-f^6 + a^6^ - 0^/3^ - a^p'^-f'^ + 4a^(3^6^ 

+ a^-f^6^ + 2aP^6 - a(5^6^ - [i^5^ + p^-f^6^) + siS2{2ap^-f5 - 2a(3j^6 - 2pS^^ 

+ 2/37^(5^) + sissi-a^p + 4a^/35 + 0^/3^ + a^p-/^ - 60^/35^ - 4a^P^6 + 40^/35^ 

+ 50/33^2 _ 3c^/3^2^2 _ ^^j4 _ 2^3 ^3 ^ 2/3725^) + slia^ - 3a^6 - a^^ - 2aS^ 

+ 80^52 ^ 2a' PH + 50^25 _ ^353 + a2^2^2 _ ^2^2^2 ^ ^2^4 _ 4^2^2^2 _ 2„^4^ 

+ 07^53 - /327252 ^ ^4^2) ^ S2S3(aS " 4a^7'^ - a^P^l - a^l^ + 6037^2 ^ 4^,2^3^ 

- 40^7^3 + 30/327(^2 _ 5c^^3^2 ^ ^^^4 _ 2/32^^3 ^ 27353) + 52^^3/32,5 - q3^2^ 

- 3a2/32(52 + 3a2^2^2 ^ 3„^2^3 _ 3^^2^3 _ ^2^4 ^ ^2^4^^ ^ ^^^^ (si(q3/3 - a^ 135 

- a(3' + aPj'^) + S2(-a^7 + 027,5 - a/327 + 073) + S3(q2/32 - 0272 - Q/32,5 + 072,5)) 
+ m^h+ (-a2/32 + 0272)] I /{2m^Lua'^{a - 6)) 
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